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We investigate simple one-dimensional driven diffusive systems with open boundaries. We are 
interested in the average on-site residence time defined as the time a particle spends on a given site 
before moving on to the next site. Using mean-field theory, we obtain an analytical expression for the 
on-site residence times. By comparing the analytic predictions with numerics, we demonstrate that 
the mean-field significantly underestimates the residence time due to the neglect of time correlations 
in the local density of particles. The temporal correlations are particularly long-lived near the 
average shock position, where the density changes abruptly from low to high. By using Domain wall 
theory (DWT), we obtain highly accurate estimates of the residence time for different boundary 
conditions. We apply our analytical approach to residence times in a totally asymmetric exclusion 
process (TASEP), TASEP coupled to Langmuir kinetics (TASEP -|- LK), and TASEP coupled to 
mutually interactive LK (TASEP + MILK). The high accuracy of our predictions is verified by 
comparing these with detailed Monte Carlo simulations. 
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I. INTRODUCTION 

Driven diffusive systems have been intensively stud¬ 
ied in the last three decades. These systems, being in¬ 
herently nonequilibrium, exhibit very rich dynamical as 
well as steady state behavior mn- A considerable the¬ 
oretical effort has been put into modeling such systems. 
One particular model that has similar appeal to the Ising 
model for equilibrium physics is the totally asymmetric 
simple exclusion process (TASEP) in ID pff^ [5111^. As 
shown in Fig. in this model a single species of particles 
performs unidirectional hopping on a ID lattice. Due to 
hard-core repulsions an empty site can be occupied by a 
single particle only. 

TASEP and its variants have been used to model intra¬ 
cellular transport m, motion of Ribonucleic polymerase 
(RNA polymerase) on the DNA (Deoxyribonucleic acid) 
during transcription and of Ribosomes on the RNA dur¬ 
ing translation [IMS]. Advances in experimental meth¬ 
ods to track macromolecules moving on ID tracks inside a 
cell |2ni - [27| have enabled a potentially direct link between 
theoretical work on the TASEP model and experiments. 

In this study we focus on the on-site residence time of 
system particles, defined as the average amount of time a 
particle spends at a particular site. Our study is directly 
relevant to the residence time of motor proteins moving 
on biofilaments inside a cell, as one can experimentally 
measure the time that an attached motor spends on the 
filament by labeling motor proteins with fiuorophores. 
The study of on-site residence times is closely related to 
the studies on dynamics of a tagged particle in an ex¬ 
clusion process |28l-[30|. In these studies the focus is on 
the distribution of a tagged particle’s position performing 
exclusion process on a ID infinite lattice. The position 
of the tagged particle in time is directly related to the 
on-site residence time of the tagged particle. In this ar¬ 
ticle we especially focus on the on-site residence times 


of a TASEP system that exhibits a first-order-like phase 
transition. This transition, which occurs for a certain set 
of boundary conditions, comprises of two phases, low and 
high density (LD and HD), which coexist on a lattice sep¬ 
arated by a shock-interface lilisiEi[ig. The interface can 
be delocalized over the entire lattice |8m33| . However, 
in a modified version of TASEP, in which the number 
of particles is not conserved on the lattice due to ad¬ 
sorption (desorption) from (to) an infinite reservoir, the 
shock interface is localized I9j[ini. This model is known 
as the TASEP coupled to Langmuir kinetics (LK) which 
we refer to as the TASEP + LK model (see Fig.[^. A fur¬ 
ther modification in which TASEP is coupled to the mu¬ 
tually interactive Langmuir kinetics (TASEP + MILK) 
can lead to even stronger localization than the TASEP 
+ LK model ElESj. Despite the interaction with the 
bulk, the ID system exhibits phase transition behavior 
under the mesoscopic scaling of kinetic rates (attach¬ 
ment/detachment rates from the lattice). By making 
the cumulative kinetic rates independent of the lattice 
size, boundaries can compete with the bulk even in the 
thermodynamic limit pun). This mesoscopic scaling of 
the kinetic rates is potentially applicable to molecular 
motors performing directed motion along ID molecular 
tracks inside cells. Typically, kinetic rates are such that 
the motors move on a significant fraction of the track 
before undergoing detachment [35]. This allows for the 
bulk dynamics to compete with the boundary, potentially 
giving rise to unusual nonequilibrium stationary states. 

The on-site residence time of a particle performing 
TASEP, when summed up over the entire lattice of a 
given size yields the total residence time of a typical par¬ 
ticle between the injection event at one boundary to the 
extraction event at the other boundary. By relating the 
residence time to the size of the lattice, one can obtain 
transport coefficients analogous to mobility and diffusion 
constant. In the presence of Langmuir kinetics, the total 
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residence time will not be equal to the sum of the on-site 
residence times as the possibility of a particle detaching 
before reaching the other end of the lattice also needs 
to be taken into account. In this study we focus on the 
on-site residence time of TASEP as well as TASEP+LK 
systems. We first present a mean-field description of the 
on-site residence time. Mean-field theory predicts that 
the residence time depends only on the steady-state den¬ 
sity of particles. However, we show that the mean-field 
estimate significantly underestimates the residence time 
of systems in LD-HD phase (the two-phase coexistence 
state mentioned previously) due to neglect of time corre¬ 
lations in the density. The time correlations, irrelevant 
in calculation of the steady-state density profile, become 
important in calculation of the residence time. We show 
that the time correlations are especially significant near 
the average location of the shock interface. We present a 
non mean-field theory of calculating the on-site residence 
time. As an input to our theory, we need the spatial dis¬ 
tribution of the shock interface in the steady state which 
we obtain using the domain wall theory [SUES]. We 
present analytical estimates of the residence time and 
demonstrate, comparing with numerics, the high accu¬ 
racy of our theory. 

The article is organized as follows. In Sec. [IT] we 
present the model considered in this study. We briefly 
summarize the findings of previous studies on TASEP as 
well as TASEP + LK which are relevant to our study. 
In Sec. |III| we present mean-field theory for calculation 
of the on-site residence time. We neglect any time corre¬ 
lations and demonstrate that the on-site residence time 
is only dependent on the steady state density profile. In 
Sec. |IV[ we present non mean-field theory for calculating 
residence time. The derivation of the non mean-field the¬ 
ory makes use of two assumptions on the movement of the 
shock interface in LD-HD phase. These assumptions are 
supported by investigating the flag model for the move¬ 
ment of the shock in Sec. |V] and extending this model 
to TASEP+LK. In Sec. |VI| we present our theoretical es¬ 
timates for the residence time and compare them with 
numerics. We apply our theoretical anlysis to TASEP 
(with and without LK). In Sec. VH we present a general 
approach to residence time by introducing a probe parti¬ 
cle and apply our findings to the TASEP + MILK case. 
In Sec.|VTTTj we present a transition from our predicition 
of residence times to the mean-field theory prediction by 
breaking one of the assumptions on the movement of the 
shock interface. Finally in Sec. |IX| we summarize our 
findings and present a brief outlook. 


II. THE MODEL 

The TASEP model is described as follows. The system 
comprises of a ID lattice of length L and N sites. Par¬ 
ticles are injected into the lattice from the leftmost site 
with a rate of a and detach from the rightmost site with 
a rate of The totally asymmetric motion of a par¬ 


ticle corresponds to its motion exclusively to the right. 
The conventional parameter p for the hopping rate on the 
lattice is set to one here. Every particle is subjected to 
hard-core repulsion. Each site is thus either unoccupied, 
or occupied by a single particle. 

The TASEP model coupled to Langmuir kinetics is 
an extension of this model, in which each particle has a 
probability of wd of detaching during a single timestep, 
and each empty site has a probability of wa of a particle 
attaching to it from the bulk reservoir during a single 
timestep (See Fig.[^. We assume that the bulk reservoir 
has an infinite capacity to act as both source and sink 
of particles. We see that the original TASEP model is 
recovered if we set uia = <^d = 0- 



Figure 1: Schematic of the TASEP model with Langmuir ki¬ 
netics. Particles are injected at the first site with rate a and 
extracted at the last site with rate /?. The motion of par¬ 
ticles is totally asymmetric as they move exclusively to the 
right with a unitary rate given that the site is empty. Hop¬ 
ping over an occupied site is not allowed. Pure TASEP is 
comprised of injection, extraction, and asymmetric hopping 
of particles interacting via hard core repulsion. In the TASEP 
+ LK model, Langmuir kinetics are additionally included by 
allowing particles to detach from an occupied site with rate 
ujd, and attach to an unoccupied site with rate uja- 

The sites are indexed from z = 0 at the injecting 
boundary to z = at the extracting boundary. The 
total length of the lattice is L and thus each site has a 
spatial extent oi e = L/N. Under the assumption that in 
steady-state spatial correlations can be ignored, one can 
derive the following mean-field equation for the average 
density p{x) of particles [51110]: 


+ {2p - l)dxp + Da(1 - p)- ^dP = 0. (1) 

In the equation above, the continuous spatial parame¬ 
ter X runs from 0 to L and corresponds to the location 
of site z through x = ie. The parameters VIa and fir) 
are referred to as the total attachment and detachment 
rates, respectively and are defined as = loaN, and 
0,0 = (jJdN. Mesoscopic scaling of the total kinetic rates 
corresponds to fixing Qa and Qo for a given L. Such 
scaling ensures that N can be varied without modifying 
the total kinetic rates from the lattice. It is only under 
this scaling that one obtains boundary driven phase tran¬ 
sition [9| by letting N —>■ oo such that e —>■ 0, reducing 
eq. 0 to a first-order differential equation. In the ther¬ 
modynamic limit of A —>■ oo, the first-order differential 
equation is overdetermined as there are two boundary 
conditions given by p(0) = a and p{L) = 1-/3. Setting 
we have three solutions to the resulting 
differential equation: 
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Pld = a + ^x (2) 

PHD = 1 — /3 — + rix (3) 

PMC =2' (4) 

The solution pld is referred to as the low density (LD) 
phase. The solution phd, as the high density (HD) phase, 
and pmc as the maximnm cnrrent (MC) phase. In this 
study we focus on shock-type profiles, which one obtaines 
when D < l — fi — a and a,fi< An example of a shock- 
type density profile is shown in fignre and |^. The 
shock-type profile is a two-phase coexistence scenario in 
which the density makes an abrnpt transition between 
Pld and phd- The mean location of the shock is de¬ 
termined by finding the point on the lattice where the 
current p(a;)(l — p{x)) of the low density phase matches 
that of the high density phase [9] . 

The original TASEP model, described in Ref. [37], is re¬ 
covered by setting D = 0. An interesting modification 
to the original TASEP + LK model has been recently 
proposed in Ref. [34]. The authors considered modified 
LK such that the local attachment and detachment rates 
of a particle depend on the the state of the neighboring 
sites. This model is referred to as the TASEP + MILK 
model and exhibits qnalitatively different phase behavior 
from the TASEP + LK model. However, one also obtains 
two-phase coexistence, as in TASEP + LK. We defer the 
discnssion of TASEP + MILK model to section |VH| 

The two phase-coexistence in TASEP + LK is strictly 
present only in the thermodynamic limit of TV —> oo along 
with the mesoscopic scaling of the kinetic rates. In pure 
TASEP, however, shock interface can exist for a finite 
N for appropriately chosen bonndary conditions. For 
any finite N, the shock interface is delocalized around its 
mean-location [3S]. Snch delocalization, in steady state, 
resnlts in a smoothing of density profile from low den¬ 
sity to high density phase. For finite N, it is interesting 
to pose the following question: if a particle is injected 
from the leftmost site of the lattice, how long does this 
particle typically stay on the lattice? And what length 
of time will it typically spend on each site? In particu¬ 


lar, would the dynamics of a delocalized shock interface 
have important implications for the residence times? Or 
can one ignore the dynamics entirely and obtain the res¬ 
idence times in a mean-field fashion? To address these 
questions, we first consider the mean-field approach be¬ 
low. 


III. MEAN-FIELD APPROACH TO RESIDENCE 

TIME 

The residence time on a lattice with N sites, as defined 
in Ref. [33] is the average time a particle entering the 
lattice spends on the lattice before detaching from the 
lattice. In this article we additionally make use of the on¬ 
site residence time: the average length of time a particle 
attaching on a particular site spends on that site before 
making a hop or detaching. For the TASEP+LK model, 
we only consider the residence times of particles attaching 
at the a boundary. 

We first consider a nnidirectional random walker on a 
one-dimensional lattice, located on a site with discrete 
index i. This index rnns from i = Q to i = N. At each 
timestep it has a probability of Pi+i to detach from site i 
(which conld either be due to hopping to the next site, or 
due to a detachment event), and a probability of 1 — 
to stay at site i. One can calcnlate the total time that a 
particle spends on site i before detaching from that site 
as follows: 

OO -| 

^n(l-pi+i)"-i(pi+i) = -. (5) 

In this snmmation, each waiting time n is weighted by 
the probability of the particle to detach after exactly n 
timesteps. 

Eq. (§ assnmes that the detachment probability is a 
Poisson distribution. An explicit expression for the aver¬ 
age on-site residence time of a particle on site i, that 
takes the system dynamics into account can be written 
as 


Ti = 


(1 — ujd )(1 — W_|_i(t)) -|- ojd + 2W_|_i(t)(l — a;£))[(l — W+i(t -|- 1))(1 — ujd ) + ujd ] 

+ 3W-i-i(t)Ai+i(f -j- 1)(1 — a;£))^[(l — 2))(1 — ujd) + wd] 


j-i 


1 — Ni+i{t){l — lod) + ^^(j + 1)(1 ~ Ai+i(t -|- j)(l — a;D))(l — LOj^y iVi+i(t -\- k)\ . 
\ j=i fc=o / 


( 6 ) 


Here, JVi(t) denotes the occnpation of the site i at time t and can be either 0 or 1. The symbol t, used throughout 
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the article, denotes discrete time and advances in units 
of 1. The average is to be understood as an average over 
several realizations of the same system in steady state. 
In the rest of this article we refer to the quantity between 
the brackets in eq. ([^ as i.e. we write 

ri = {Ti{t)}. (7) 

In order to obtain a closed form expression for the 
on-site residence time, the time correlations of the form 

n 

(n + k)) are needed for all n G N. Under the 
k=0 

assumption that the time correlations can be neglected, 
i.e., 

n n 

(IT + *)) = n + fc)) = 

fc=o fc=o 

the above expression simplifies considerably and can be 
expressed explicity in terms of the steady state density 
Pi+i as: 


This reduces to o 1 — pi+i for the pure TASEP 
model, which is expected in the mean-field description 
|39| . Similarly the boundary condition on the rightmost 
site N is given by 

(1 - ujd)I3 + lud ' 

For the pure TASEP model, the total residence time 
7?tasep can simply be written as 

N 

7?tasep = ^ ri . (10) 

i=l 

For the TASEP+LK model however, the possibility of 
a particle detaching before reaching the other end also 
needs to be taken into account. The modified equation 
for R takes the following form: 


ri = 1- pi+i[l - ujd) + 2/9i+i(l - wd)[1 - Pi+i{l - wd)] 
+ 3[/3i+i(l — — iod )] + 

OO 

= XI - P*+l(l - ujd)) 

n—1 


1 - /Oi+i(l - ujd) ' 

The on-site residence time in eq. (|^ is the same as in 
eq. (|^ under the identification of pi+i GG 1—pi+i(l— 


N i 

7?tasep = ri + E^n fk-l,k- (11) 

i=2 k=2 


For this modification, we consider the ensemble of steady 
state configurations that have a particle occupying the 
site {i — 1). The factor in eq. (Ill, defined for 

2 < i < N, denotes the probability that a particle in this 
ensemble hops from site {i — 1) to the next site i (does 
not undergo detachment). This probability is given by 


= ( (1 - ujd)[1 - A^,(t)] + (1 - uJD)^N,{t)[l - N,{t + 1)] 

-f (1 — ujd)^N i{t)Ni{t + 1)[1 — Ni{t 2)] + ...^ 

/ OO f -1 \ 

= ( (1 - ujd){ 1 - N,{t)) + E(1 - wb)^-'(1 - N,{t + j)) J] N,{t + k)). 


1=1 


k=0 


( 12 ) 


Applying the same mean-field approximation as in the 
derivation of eq. (|^, the equation for can be writ¬ 

ten as 


OO 

= Epr'(i-^iir(i-pO 

n—1 

_ ^ __ 

1-Pj(1-w_d)' 


(13) 

(14) 


To calculate the total residence time for the TASEP+LK 
model, we multiply the residence time of each site by the 
probability that a particle starting at the a boundary 
reaches this site (and thus does not detach before that). 
We thus write 
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^TASEP+LK = ■'’1 + 1^1 — 
N 


UID 


r’2 + 1 - 


LUd 


1 — ^2(1 — ^ d ) 

n , 


n=2 fe=2 




1 - P2{1 - LOd) 
N 1 


1 - 


UJd 


1 - P3(l - Wd) 


rs + ••• 


= 'Tl + ^ r„ (1 - LOoTk)- 


(15) 


n=2 k=2 


We see that eq. (151 reduces to eq. (101 if we set ujd 

0 . 


A. Breakdown of mean-field in LD-HD phase 


Eqns ([^ , (10) and (151 are valid if the time correlations 
in local density are negligible. If, however, there are time 
correlations present in the local density, the mean-field 
expression underestimates the on-site residence times. 
Monte Carlo simulations revealed significant time cor¬ 
relations for systems in LD-HD phase, see fig. |^e & f. 
The time correlations C'i(At) are calcnlated by 


+ At) - p,)) 


(16) 


where each Monte Carlo timestep in the random sequen¬ 
tial update (for a description of this update procedure, 
see e.g. sni) corresponds to a time interval of At = 1. 
The time correlations are rather long-lived in compari¬ 
son to the mean-field estimate of residence time, espe¬ 
cially in the neighborhood of the average shock location. 
We therefore expect eq. (|^ to give incorrect values for 
the on-site residence times for these profiles. In order to 
obtain accnrate on-site residence time, a non-mean-field 
expression for for these profiles is thus needed. 



Figure 2: Left panels correspond to pure TASEP and the right 
panels to TASEP + LK. (a) and (b): Steady state density 
profile, (c) and (d): On-site residence times from simulations 
(dashed lines) together with the mean-field prediction (solid 
lines), (e) and (f): Time correlation of local density at dif¬ 
ferent time steps (see legend). Model parameters for TASEP 
+ LK: a = P = 0.1, fl = 0.1 and N = 100. TASEP model 
parameters: a = P = 0.1 and N = 100. 


IV. NON MEAN-FIELD APPROACH TO 
RESIDENCE TIME 

A shock-type profile is characterized by the existence 
of a low-density phase at the a bonndary and a high- 
density phase at the /3 bonndary. The sharp transi¬ 
tion from one density phase to the other is considered 
as the shock interface. Two plots of the density pro¬ 
files of snch LD-HD phases are shown in Fig. In the 
TASEP model, the shock is completely delocalized. In 
the TASEP+LK model, the interface is localized only in 
the limit of A —)■ 00 . For finite A, the interface per¬ 
forms a random walk in a confining potential [38]. As 
shown above, time correlations in the local density are 
particularly significant near the average shock location 
persisting on time scales longer than the mean-field es¬ 
timate. In order to derive an expression for the on-site 
residence times in presence of temporal correlations, we 
make the following assumptions about our system: 


1. Each system site can at each timestep be deter¬ 
mined to be either in the low density phase or in 
the high density phase. 

2. A particle waiting to hop on site i stays in the same 
density phase until it has hopped. 

In the following section. Sec. |V A] we investigate the va¬ 
lidity of these assnmptions. 

We start by rewriting eq. (§ as 

Ti(t) = [$i(f) -b (1 - ^>*(0)] T-i(t), (17) 

where we define $i(t) as a time-dependent phase con¬ 
stant. $i(t) takes on value of 1 if site i belongs to the 
low-density phase and 0 if site i belongs to the high- 
density phase. 

We now take the average of both sides of this eqnation 
in the following way: 

{n{t)) = {^i{t)n{t)) + ((1 - $,(t))Ti(t)). (18) 
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Each of the two terms on the right hand side corresponds 
to an ensemble average. Using assumption 1, the fraction 
of systems in the ensemble that are in the low-density 
phase on site i is given by ($i(t)). Similarly, the fraction 
of systems in the ensemble that are in the high-density 
phase on site i is given by (1— <f>i(t)). However, since ^i{t) 
is 1 in the low-density phase and 0 in the high-density 
phase, the term ($i(t)Ti(t)) in eq. (18) gets a contribution 


only from those systems in the ensemble which are in 
low-density phase. Similary, ((1 — $i(<))Ti(t)) gets finite 
contribution only from the systems in high-density phase 
in the ensemble. We can thus write 

{n{t)) = initMiMt)) + (19) 

where and (...)f° denote an average over the 

systems in LD phase and the systems in HD phase at 
site i at time t respectively. 

We now have a weighted sum of the average residence 
time in the low-density phase and the average residence 
time in the high-density phase. In a recent study [3S] 
it has been shown that the mean-field prediction for the 
residence times is accurate for a purely HD system as well 
as a purely LD system. We now make use of assumption 
2 above, which allows us to identify with the 

on-site residence time on site f of a pure LD system, and 
with the on-site residence time on site i of a pure 
HD system. We thus arrive at the following expression 
for the on-site residence time for an LD-HD system: 


A,; 


1- A,; 


Ti = 


1 — PLD,i(l — Wd) 1 — pHD.i(l — 


( 20 ) 


Here, we have identified ($i) = A^. 

As our parameter Ai is equal to the fraction of time that 
site i is part of the low-density phase, and 1—Ai is equal to 
the fraction of time site i is in the high density phase, we 
can write the following expression for the average density 
profile: 

Pi = AiPLD.j + (1 — Ai)pHD,i. (21) 

From this, we can obtain an expression for Ap 


A,; = 


Pi — PHD,i 
PLD,i — PHD.i 


( 22 ) 


This equation, together with eq. (201, gives us a com¬ 


pletely analytical result for the on-site residence times in 
terms of pld.i, Phd.i and pi. 

Eq. 

fP in LD-HD phase. For TASEP+LK how¬ 


(101 for the total residence time R still holds for 
pure pfe" 


ever, eq. (151 is no longer valid as eq. (13) ignores time 


correlations. A non mean-field expression for fi-i^i for 
systems in LD-HD phase can be derived using the same 
approach as in the derivation of eq. (20). The result is: 


fi-i,i —1 — 


uid\ 


1 — PLD.i(l — Wd) 
LUpil — Xj) 

1 — PHD,i(l — W_d) 


(23) 


The total residence time R on the lattice is then given 
by 


N 


Rtasep+lk = ri + ^ ^ ~ T 


UlD^k 


n=2 k=2 - --PLD.fc(l-a;D) 

w d(1 — A/c) 


1 — PHD,fc(l — Wd) 


N n 

= J'l + ^ (1 - ujork)- 

n—2 k—2 


(24) 


This equation is identical to eq. (15), but with now 
given by eq. (20). 


A. 


V. DOMAIN WALL DYNAMICS 

Flag theory for TASEP+LK with random 
sequential update 


In the previous section, we saw there are two conditions 
that need to be satisfied in order for eq. (20) to hold: (1) 


a location can be assigned to the shock at each timestep 
and (2) the shock moves in such a way that a particle 
crosses the shock only when it makes a hop. In order 
to investigate if these conditions are indeed satisfied, a 
description of the shock between the low density phase 
and the high density phase is thus needed. 

For the pure TASEP model, a microscopic description 
of the shock location has been derived in Ref. [IT]. This 
was done by introducing a second class particle. This 
is a particle that hops like the system particles if it has 
a hole to its right, but if a system particle to its left 
tries to hop, the particles switch places. In Ref. [H], it 
was proven that this second class particle indeed tracks 
the location of the shock. This implies that for the pure 
TASEP model, our two conditions are indeed satisfied: 
all sites can be determined to be either in the high or the 
low density phase, and a particle makes the transition 
between phases only when it hops onto the site where 
the second-class particle is located. 

We would like to extend the idea of this shock marker 
to a TASEP+LK system. A problem for doing this is 
that this second class particle would have to be subject 
to langmuir kinetics, allowing it to detach, and thus re¬ 
moving the shock marker from the system. 

A solution to this is to consider a flag as described in 
Ref. |12] instead of the second class particle. This flag 
was introduced in Ref. |12] to track the location of the 
shock in pure TASEP with parallel update, and obeys 
the following motion rules: 

1. The flag is first placed at the leftmost particle of 
the leftmost cluster 

2. If the particle with the flag on it moves, the flag 
moves with it 
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i i+1 i+1+m i+mt-2 

Figure 3: A graphic representation of flag theory, (a) A flag is 
placed on a system particle to keep track of the shock location. 
Its movement rules are such that the flag is always located 
on the leftmost particle of a cluster, (b) and (c) TASEP 
+ LK: If two clusters merge forming a single cluster due to 
an attachment event, the flag which was previously on the 
right cluster, moves to the leftmost particle of the new merged 
cluster. In the event of detachment of the particle carrying 
the flag, the flag is transferred to the leftmost particle of the 
cluster immediately to the right of the detaching particle. 

3. If a particle is blocked by the particle carrying the 
flag, the flag is transferred to this blocked particle. 

If we apply these motion rules to a TASEP system with 
random sequential update, we see that the motion of the 
flag is identical to the motion of the second class particle. 

A note on the application of rule 1 in our model: due 
to clusters spontaneously forming in the low density re¬ 
gion, the flag is not necessarily placed close to the shock 
interface using this rule. However, it ensures that the flag 
is placed on particle such that the flag has no particle to 
the left of it. After a few timesteps, the flag will move 
close to the shock interface due to rules 2 and 3 above 
and will be a good indicator of the shock position. 

For the TASEP+LK and TASEP+MILK cases, we 
have additional movement rules for the flag due to attach¬ 
ment and detachment events. We propose the following 
additions (shown graphically in Fig.|^: 

4. If the flag is located at site i, a particle attaches 
at site i — 1 and sites i — 2 — n through i — 2 are 
occupied, the flag moves to site i — 2 — n. If site 
i — 2 is unoccupied, the flag moves to site i — 1. 

5. If the flag is located at site i, the particle at site 
i detaches and sites i + 1 through i + 1 + m are 
unoccupied, the flag moves to site i + l + m+1. If 
site J -I- 1 is occupied, the flag moves to site i + 1. 

These rules are natural if an attachment event between 
the flag and a cluster ending two places to the left of 
it is viewed as several particles coming in at the same 
timestep. Similarly, a detachment event of the particle 
which carries the flag and has several holes to its right is 
viewed as several holes coming in at the same timestep. 

To test the validity of this description of the flag loca¬ 
tion, we performed simulations keeping track of the flag 
locations together with the average density profile. The 
distribution of flag locations was used together with the 



Figure 4: Average density profiles for TASEP+LK obtained 
by keeping track of the shock location using our rules for 
the flag movement. The solid red line corresponds to the 
profile predicted from the flag locations, the blue dashed line 
corresponds to the actual density profile obtained in the same 
simulation. The density profiles of figures a), b) and c) are 
taken over a relatively short timescale (600 timesteps in the 
Monte Carlo simulation), the profile of figure d) is taken over 
a long timescale (50000 timesteps). The excellent agreement 
between the predicted and numerically obtained profiles in 
(d) is a clear indication of the accurate recovery of shock 
position distribution in the steady state. Model parameters: 
a = l3 = n^ 0.1, N = 100. 

LD and HD average density solutions eqs. and (i to 
predict the average density at each site. Density profiles 
obtained this way matched very well with the actual av¬ 
erage density, for short as well as for longer run times, 
as shown in Fig. With the above modifications, flag 
theory can thus be applied to TASEP+LK as well. 

With these additions, it is now possible for a particle 
to switch phases while waiting to hop. The derivation in 
the previous section is thus no longer strictly valid. It is 
however still a good approximation if such a switch is a 
rare event. We therefore consider the fraction x of flag 
movements due to attachment and detachment events. 
For a flag located at site i, this fraction is given by 

X — {Ni-2{1 — Ni-i)uJA + NiUJuil — -^i+l)) 

= ujA{N,-2il - fV,_i)) + - N,-i))- (25) 

We see that this fraction is much smaller than 1 in the 
case that uja^ojo 1- This is indeed the case for LD-HD 
systems as long as the number of sites N is not of order 
1. To see this, we note that for the TASEP+LK model, 
we have a shock in the case that H < 1 — /3 — a. This 
implies that uja,^d < so uia,<^d can only be of 

order 1 if is of order 1. 

We conclude that the two assumptions for the deriva¬ 
tion of the previous section are indeed valid. It is possible 
to assign a location to the shock at each timestep, thus it 
is possible to determine each system site to be either in 
the low density phase or in the high density phase at each 
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timestep. As long as the number of sites N is not of or¬ 
der 1, a particle, to a good approximation only switches 
phases while performing a hop. Thus it is indeed justi¬ 
fied to write eq. (201 for the on-site residence times for 
LD-HD systems. 


B. Density profiles in LD-HD phase 

In section EYl eq. 1^ for the on-site residence times 
for an LD-HD system was derived in terms of pi^d{x), 
PHD (a;) and p{x). In order to calculate on-site residence 
times from this equation, an equation for p{x) is thus 
needed. A derivation for this has been done in Ref. [3H] 
using Domain wall theory (DWT). We briefly summarize 
their results below. 

The average density profile in this system (written in 
the continuum limit) is described by 

p{x) = phd{x)[ p{x')dx'+ plu{x) f p{x')dx'. (26) 
Jo Jx 

In this equation, p{x) describes the distribution of the 
locations of the shock. The system length L has been set 
to 1. The solution for p{x) found by (38] is given by 


p{x) = 


-N R (1-5 

o Jccn ' ( 


NuJi{x) 




)dx 


(27) 


Here, N is the number of sites on the lattice, and ujr and 
oji are the hopping rates of the shock location to the right 
and left respectively as shown in Ref. [55] 

PHD 


UJi = 


PHD — PLD 
Pld 

PHD — Pld 


The normalization constant N is chosen such that 
Jp p{x)dx = 1. The term xq can be chosen arbitrarily 
as it will be absorbed in the normalization constant N. 

In the following sections, we make use of these results 
together with our eqs. (201 and (221 to obtain residence 


time predictions in specific conditions. 


VI. RESULTS 


A. Shock-type profiles in the TASEP model 


For the pure TASEP case, we have a completely de- 



Figure 5: On-site residence times for different system lengths. 
The left panels correspond to pure TASEP with a = (3 = 0.1, 
the right panels correspond to TASEP+LK with a = P = fl = 
0.1. Figures (a) and (b) correspond to A = 50, (c) and (d) 
to A = 100, (e) and (f) to A = 200. Green dashed lines: on¬ 
site residence times obtained from simulations. Black solid 


lines: analytical on-site residence times from eq. (20l. Red 


dash-dotted lines: mean-field result for residence times. 


We see here that our distribution function is no longer a 
function of the number of sites N. This means that the 
shape of our profile is now size-independent. 

Using eq. (261 in the continuum limit, the expression for 


the average density profile p(x) becomes 
p{x) = a + x(l — 2a), 
which agrees with the result in Ref. j^. 


(29) 


Using eq. (201 to obtain the on-site residence time we 
obtain 


r(x) = 


a + x(l — 2a) 
q:(1 — a) 


(30) 


The results from Monte Carlo simulation together with 
the prediction by eq. (301 are shown in Fig. The total 
residence time on the lattice is now given by 


— -, - - ^ 

localized shock in the case that a = j3 < Our two A = A / r{x)dx = N 

solutions for the differential equation (fll that meet a Jo 


solutions for the differential equation p 
boundary condition are given by eqs. (]2 
0 = 0. 


that meet a 
and (§ with 


2q;(1 — a) 


(31) 


We use eq. (271 to calculate p{x). Since Ur = uji for 
the pure TASEP model, this expression reduces to 


We thus recover the linear dependence of the residence 
time on the system length, in accordance with the re¬ 
sults of Ref. (Ml. The results for our numerical simu- 


Ptasep{x) — 1. 


(28) 


lation compared to the results predicted by eq. (311 are 
presented in Fig. 
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B. Shock-type profiles in the TASEP+LK model 

The derivation of the analytical profile for a shock-type 
profile in the TASEP+LK model has been done explicity 
in Ref. [3S] for the case that 17^ — fli) — fl. The result 
is: 


P(^) 


1 + erf 2- 


NnA 


(1 + A)(1-A) 


(x - Xs) 


+ fix + a, 


(32) 


where A is the height of the shock, given by A = phd — 
Pld, and Xs is the average shock position. The on-site 
residence time can now be calculated using eq. (201 in 
the continuum limit: 


1 — {fix + q:)(1 — ojd) 

1 — A(a;) 

1 — {fl{x — 1) + 1 — /3)(1 — ujd) ’ 
where X{x) is given by 


(33) 


X{x) = 


1 + erf 


(V 


NQA 


(1+A)(1-A) 


^{x-x, 


1 -\- a — fl — (3 


1. (34) 


The results for our numerical simulation compared to 
the results predicted by eq. (331 are presented in Fig. ^ 
The total residence time R can be obtained numerically 
using eq. (241 with eq. (331 as arguments. The mean-field 
result would be given by evaluating eq. (15) with eq. ^ 
as arguments. The results of monte-carlo simulations are 
shown in Fig. For large N, the iV-dependence of R is 
very close to linear as can be seen in Fig. An explicit 
expression for R in the large N limit is given in appendix 
[B] eq. (B7). Since i? ^ A, it follows that one can define 
a transport coefficient analogous to mobility (a material 
property independent of system size) in the limit of large 
N for TASEP as well as TASEP + LK. 


Until now, we have focused on systems with boundary 
conditions chosen such that the steady state density pro¬ 
file is described as the LD-HD phase. In the following 
section, we consider systems in MC phase and demon¬ 
strate that our theoretical approach is equally applicable 
to these systems as well. 



N N 


Figure 6: Total residence times R for TASEP in LD-HD phase 
(a), for TASEP+LK in LD-HD phase (b) and for TASEP+LK 
in MC phase (c). Blue dash-dotted line: results from Monte 
Carlo simulation. Black line: non mean-field prediction. 
Green longitudinally-dashed line: mean-field prediction. Red 
transversally-dashed line: prediction by eq. (B71. Panel (a): 
non mean-field result calculated with eq. (101, model parame¬ 
ters a = P = 0.1. Panel (b): non mean-field result calculated 
with eq. (2^, model parameters a = P = fl = 0.1. Panel (c): 
analytical result calculated with eq. (361, model parameters 
a = P = 0.5, fl = 0.1. 


In a maximum current phase, characterized by a con¬ 
stant density in the bulk, depending on the mismatch 
between the boundary conditions and the Langmuir 
isotherm, boundary layers can exist. The width of the 
boundary layers relative to the system size approaches 
zero in the limit of A —oo. Below we assume that the 
boundary layers can be neglected in calculating the total 
residence time. 

An exact expression for the total residence time for a 
constant density can be obtained by inserting = p in 
eq. (24) which yields 


R = 


{N-fl){l-p)(l- 


1 - 


a 


Qp—pN-\-N 


N 


U(l-p(l-R)) 


(35) 


to 


It can be easily shown that for N ^ eq. (351 reduces 
R= N/fl. (36) 


C. Residence time for TASEP+LK in maximum 
current phase 

The residence time of a particle when the system is in 
MC phase has been studied in Ref. [35]. It was shown 
that the total residence time scales linearly with the 
system size. However, the study focused only on pure 
TASEP. We will now show that this linear scaling re¬ 
mains present for a TASEP+LK system as well. 


It follows that for a constant bulk density, the total resi¬ 
dence time scales linearly with the system size. The fac- 
tor e i-p is the fraction of particles which starting on the 
first site of the lattice hop all the way to the /3 boundary. 
This factor is 1 in the limit of fl/{l—p) +; 1 in which case, 
the total residence time scales asi?«A/(l — p). On the 
other hand, when 0/(1 — p) ^ 1, almost every particle 
detaches from the lattice before reaching the /3 boundary. 
In that case, the total residence time, consistent with the 
assumed mesoscopic scaling is Rw I/ojd = N/fl. 
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VII. RESIDENCE TIME OF A PROBE 
PARTICLE 

Including mutual interactions between particles in ad¬ 
dition to the hard-core repulsion gives rise to several in¬ 
teresting features. In the Katz-Lebowitz-Spohn (KLS) 
model, the hopping rates are modified depending on the 
occupancy of the next-nearest neighbor resulting in ad¬ 
ditional correlations |13]. This model gives rise to exotic 
features such as localized downwards shocks and phase 
separation into three distinct regimes [H]. In the model 
TASEP + MILK shown schematically in Fig. the at¬ 
tachment and detachment rates are modified depend¬ 
ing on the occupancy of the adjacent sites [51]. These 
additional interactions increase (decrease) the effects of 
boundaries on the phase behavior of the model. For such 
models with additional interactions, it is in principle pos¬ 
sible to calculate the on-site residence times. However, 
a more general and considerably simple description of 
residence time is possible by introducing a probe parti¬ 
cle. The probe particle interacts via hard-core repulsion 
and does not undergo detachment. We consider a system 
consisting of of a single probe particle and all the other 
particles following TASEP (with or without LK). It is 
interesting to study the residence time of a single such 
probe particle. 

A. Probe particle dynamics 

The residence time of the probe particle can be calcu¬ 
lated using the same approach as in the derivation of eq. 
^ for ri in the TASEP+LK case. The difference is now 
that the probe particle cannot detach. As the w^i-factors 
take the detachment of the particle considered into ac¬ 
count, the on-site residence time of the probe particle is 
obtained by setting ujd = 0 in eq. (|^. 

Note that the resulting expression is a general one, 
valid for a probe particle in the TASEP model with any 
type of further modified interactions. The modifications 
of the interactions will all be captured in the average 
occupations and the time correlation functions. If time 
correlations can be ignored, this expression reduces to 


recently proposed in |34| : the totally asymmetric exclu¬ 
sion process together with mutually interactive Langmuir 
kinetics. 


B. Results for a probe particle in a TASEP+MILK 
system 

A schematic of TASEP + MILK system is shown in 
Fig. [11 Every particle is subjected to the same rules as in 
TASEP. In addition, every particle can undergo attach¬ 
ment and detachment as in Langmuir kinetics with the 
modification that the kinetic rates are dependent on the 
state of neighboring sites. If a particle has one neigh¬ 
bouring particle, its detachment rate becomes yw/j. If it 
has two neighbouring particles, its detachment rate be¬ 
comes 7 ^ 0 ;£). If there are no neighbouring particles, its 
detachment rate remains ujd- Similarly, the attachment 
rate becomes Suja or 6'^uja if one or both neighboring sites 
are occupied, respectively. If there are no neighbouring 
particles, its attachment rate will remain equal to u>a- 
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Figure 7: A graphic representation of the TASEP+MILK 
model. Every particle performs ordinary TASEP in addition 
to Langmuir kinetics which depend on the ocupation states of 
neighboring sites. The detachment rate of a particle is scaled 
by a factor 7 if one of its neighboring sites is occupied. If 
both neighboring sites are occupied then the rate is scaled by 
a factor of 7 ^. Similarly, the attachment rate is scaled by S 
when only one of the neighboring sites is occupied and by 5^ 
when both are occupied. 

The differential equation describing the steady-state 
solutions density profile to first order in e is given by [31] 


1 



(37) 


^^'xP A (2p ~ ^)dxP + Ha)! + p{5 — 1)]^(1 — p) 

-Hb[1 + p(7-1)]2p = 0. (39) 


For any system in LD-HD phase subject to the two 
assumptions stated in section |TVj the on-site residence 
time is given by 


\i ^ 1 Ai 

1 — PLD 1 — PHD 


(38) 


The derivation of this equation is completely analogous 
to the derivation in section CYl 

We now apply the description of the residence times 
of a probe particle to a model with modified interactions 


This equation has not been solved in general, but solu¬ 
tions are known for certain parameter regimes |34| . 

In the TASEP+MILK model, shock-type profiles have 
been found to occur as well [34] . At the end of this section 
we apply our method to determine the residence time for 
these profiles, and compare with simulations. 

First we investigate if it is justified to use the assump¬ 
tions stated in section iTVl to the TASEP+MILK model. 
For this model, we can still use the rules for flag move¬ 
ment in a TASEP+LK stated in section |V A| This is be¬ 
cause there are no new ways for particles to hop or attach 
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Figure 8 : Average density profiles for TASEP+MILK ob¬ 
tained by keeping track of the shock location using our rules 
for the flag movement. The solid red line corresponds to the 
profile predicted from the flag locations, the blue dashed line 
corresponds to the actual density profile obtained in the same 
simulation. The density profiles of figures a) and c) are taken 
over a relatively short timescale (600 timesteps in the Monte 
Carlo simulation), the profiles of figures b) and d) are taken 
over a long timescale (50000 timesteps). Model parameters 
for a) and b): <5 = 7 = 0, a = /? = 0.15, Q = 0.1, N = 100. 
Model parameters for c) and d): 5 = l + i/j and — 1 —ip with 
Ip = 0.2, a = 0.05, P = 0.15, Q. = 0.3, N = 100. As in Fig.|4] 
the excellent agreement between the predicted and numeri¬ 
cally obtained profiles in b) and d) is again a clear indication 
of the accurate recovery of shock position distribution in the 
steady state. 


to sites in the TASEP+MILK model; the only modifica¬ 
tion is in the probability of these events. We can calculate 
the fraction of flag movements due to attachment and de¬ 
tachment events using the same approach as in section 
|V A[ the result is 


X = (40) 


We see that this fraction is much smaller than 1 as long 
as ^ 1 and Jwd +; 1. To test the validity of 

this description of the flag location in TASEP+MILK, 
we again performed simulations keeping track of the flag 
locations together with the average density profile. The 
approach is the same as for the plots in Fig. Density 
profiles obtained this way matched very well with the 
actual average densities, for short as well as for longer 
run times, as shown in Fig. 

The solutions for the density profile are only known 
in specific regimes. We therefore compute the residence 
times in each of these regimes separately. We first con¬ 
sider the the parameter regime S = In this case 
both attachment and detachment rates are symmetrically 
modified such that both the attachment and detachment 
rates are either scaled up or scaled down. The solution to 
the eq. (391 derived in Ref. 1341 is reproduced in the Ap¬ 


pendix. Using eqs. (A1) and (A2| together with eqs. (27) 


and (26) we obtain the density profile for a given sys¬ 
tem size N. We use numerical integration to obtain the 
on-site residence time from eqs. (20) and (22) which is 
shown in Fig. 

Another interesting parameter regime is S = l + ip,"f = 
l — ip. In this case mutual interactions modify the kinetic 
rates in an asymmetric fashion such that one is scaled up 
and the other is scaled down. Assuming that ip 1 
one can obtain closed form expression for the solution to 
eq. (39) as shown in Ref. [34] . The expressions are repro¬ 
duced in the Appendix (eqs. (A3) and (A4)). Performing 
similar calculations as in the case of symmetrically mod¬ 
ified kinetic rates, we obtain on-site residence time as 
shown in Fig. 



Figure 9: Residence time of a probe particle in presence of 
system particles undergoing TASEP + MILK, (a) Symmetri¬ 
cally modified kinetic rates S = 'y = 0,a = P = 0.15, Q — 0.1, 
N = 100. (b) Asymmetrically modified kinetic rates <5 = l + ip 
and 7 = 1 — '0 with ip = 0.2, a = 0.05, P = 0.15, fl = 0.3, 
N = 100. Green dashed lines: residence time of the probe 
particle from simulations. Thick solid lines: Our theoretical 
prediction. Red dash-dotted lines: Mean-field prediction of 
the on-site residence time. 

As can be seen in Fig. our theoretical predictions 
are in excellent agreement with the numerically obtained 
on-site residence time. We believe that using the concept 
of probe particle, our theoretical approach can describe 
the residence time in presence of interactions other than 
the MILK. 


VIII. TRANSITION TO MEAN-FIELD THEORY 
IN LD-HD PHASE 

The approach for calculating on-site residence times 
presented in this article is valid for any two-phase system 
that obeys the two conditions described in section jlVl To 
explore a system for which this approach does not work, 
we introduce a ghost particle to our system which we give 
the following properties: 

• The particle is first placed on the first site. 

• The particle moves to the next site if this site is 
unoccupied. 

• Other system particles are not influenced by this 
particle. The system particles can thus hop onto a 
site if it is occupied by the ghost particle, as long 
as it is not occupied by another system particle. 
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Figure 10: A graphic representation of the ghost particle dy¬ 
namics. The ghost particle is first placed on the first site. 
When in active state, the ghost particle moves exclusively to 
the right only if the next site is unoccupied. In the inac¬ 
tive state, the ghost particle always remains at the same site. 
The probability of being in active state is denoted by /. On 
reaching site N, the ghost particle detaches with a rate of /3. 
Effectively, the ghost particle follows TASEP in active state. 
Other system particles do not interact with the ghost particle. 
A site can be occupied by a system particle together with the 
ghost particle. 



Eigure 11: Plot of the on-site residence time of the ghost 
particle for different values of ghost frequencies /. The dashed 
line and the dotted line indicate the on-site residence times 


predicted by eq. (20 ) and by mean-held theory respectively. 


The dashed line coincides with the residence time of the ghost 
particle for / = 1. 


IX. DISCUSSION AND CONCLUSION 


• When the ghost particle reaches site N, it has a 
probability of (3 of detaching. 


Furthermore, we will give this particle a probability to 
pause, a concept introduced in Ref. |45| . When a particle 
is paused, it will not be able to hop. As long as the 
particle is not paused, it will move according to the rules 
above. We give the particle a probability of / to be in the 
active state on a single timestep and thus a probability 
of 1 — / to be in the paused state. We keep track of the 
on-site residence times of this particle over all timesteps 
that it is in the active phase. 


Due to these movement rules for our ghost particle, it 
is now possible for the flag described in section V A to 
pass under the second-class particle while it is waiting to 
hop. The second assumption of our approach, that a par¬ 
ticle only switches between phases while making a hop, 
will thus not be satisfied if / < 1. We thus expect the 
prediction of our theory for the crossing time of the ghost 
particle to be inaccurate in this case. If the probability 
to be in the active state / is low enough, the average 
time intervals between active states of the paused par¬ 
ticle will be larger than the time interval over which we 
observe time correlations. We then expect for approxi¬ 
mation ([^ to be valid again, and thus for the mean-field 
expression for the on-site residence times to be accurate. 
In the intermediate regime, we expect both approaches to 
be inaccurate, and we expect the on-site residence times 
to be somewhere between the mean-field value and the 


value predicted by eq. (201. 


We ran simulations for the on-site residence times of 
this ghost particle for various values of /, for two LD- 
HD profiles. The results are shown in Fig. [m As can 
be seen in this figure, the on-site residence times for the 
ghost particle indeed make a gradual transition from the 
prediction of eq. (201 to the mean-field prediction for de¬ 
creasing /. 


In this paper, we study the average residence time of 
a particle on a given site in a ID driven diffusive sys¬ 
tem with open boundaries. The particle performs to¬ 
tally asymmetric motion with hard core repulsion. In a 
mean-field scenario, ignorning spatial and temporal cor¬ 
relations, we obtain the on-site residence time as inverse 
of the average hopping rate {1 — pi{l — ojd)), where pi is 
the steady-state density of particles at site i and ujd is the 
single-site detachment probability. Using Monte Carlo 
simulations we show that the mean-field significantly un- 
derstimates the residence time in the neighborhood of 
shock interface. The shock interface is characterized by 
a steep transition in the density profile of particles in the 
steady state. In the thermodynamic limit, the transition 
is infinitely steep, reminiscent of a first order phase tran¬ 
sition with two-phase coexistence. The underlying reason 
for failure of the mean-field estimate for residence time 
is the neglect of time-correlations in the local density of 
particles. The temporal correlations in the density are es¬ 
pecially long lived near the shock interface with the local 
correlation time much larger than the mean-field on-site 
residence time. 

The shock interface separating the low and high den¬ 
sity phases is not localized on a given site. In pure 
TASEP, the interface is delocalized over the entire lat¬ 
tice for any system size. On adding LK, the interface 
becomes localized within a region that scales with \/N 
where N is the system size. For any finite size lattice, the 
interface can be considered as performing random walk in 
a confined potential. The mean-field residence time does 
not take into account the fact that the density at a local 
site fluctuates on time-scale longer than the mean-field 
estimate. One can obtain the steady state density profile 
at a given site by averaging over time scales longer than 
the correlation time. However, the presence of long-lived 
correlations leads to an underestimation of the residence 
time using mean-field. We provide a non-mean-fleld ex¬ 
pression for the on-site residence time which takes the 
movement of the shock interface into account. Our de¬ 
scription requires calculation of a single parameter, re- 
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ferred to as the site-dependent, average phase constant. 
Our derivation of a non mean-field equation for the on¬ 
site residence time relies on two assumptions: ( 1 ) that 
each site can at each timestep be determined to be either 
in the low or the high density phase, ( 2 ) that a particle 
crossing the lattice switches phases only while making a 
hop. These two assumptions are validated by consider¬ 
ing Flag theory for the movement of the shock. In previ¬ 
ous studies. Flag theory has been successfully applied to 
pure TASEP to track the shock location. Here we extend 
the Flag theory to TASEP + LK and determine the flag 
movement rules. These flag movement rules validate the 
assumptions of our non mean-field derivation of on-site 
residence times. 

We obtain the average phase constant making use of 
previous results for the average density profile from do¬ 
main wall dynamics and compare the analytical residence 
time with the numerical simulations. We show that our 
analytical predictions are highly accurate in describing 
the residence time in TASEP as well as TASEP+LK. 

The total residence time on a pure TASEP lattice can 
be obtained as the sum of the on-site residence times. 
In TASEP+LK, the possibility of a particle detaching 
before reaching the /3 boundary also needs to be taken 
into account. Doing this yields a weighted sum of the 
on-site residence times as an expression for the total res¬ 
idence time. We demonstrate that the total residence 
time asymptotically scales linearly with the system size. 

Our approach to calculating residence time can be ex¬ 
tended to TASEP systems with further modified dynam¬ 
ics by considering a test particle. This is a particle that 
does not detach and will thus cross the entire lattice. 
The on-site residence time for this particle is obtainable 
provided phDi Pld and p(^x) are known. We run sim¬ 
ulations for the on-site residence time of a test particle 
in the TASEP+MILK system and And good agreement 
with our predictions. 

In order to demonstrate how mean-field theory can be 
recovered by violating the second assumption of our ap¬ 
proach, we consider a ghost particle on the lattice which 
performs pure TASEP. The ghost particle has a hopping 
rate that is lower than that of the system particles. The 
ghost particle, somewhat similar to a second class parti¬ 
cle, is invisible to other particles but is itself subjected to 
hard core repulsion. We demonstrate that when the hop¬ 
ping rate of the ghost particle is sufficiently low, it can 
experience phase-change without performing a hop. The 
phase change occurs due to the shock interface moving 
past the ghost particle while it is paused. In the case of 
vanishing hopping rate, we find that the residence time 
of the ghost particle is accurately captured by the mean- 
field estimate. 

Our focus in this study is the residence time on a lat¬ 
tice of given size with given boundary conditions. It will 
be interesting to explore the same in a network setting, 
viewed as collection of such ID systems. The bound¬ 
ary conditions in a network will be specified only on 
the outermost boundary with all the internal nodes ex¬ 


hibiting fluctuations in the injection and extraction rates. 
Such a study could be linked to existing work considering 
TASEP on networks [IS] and TASEP+LK on networks 
m- Another interesting study would be to apply the 
concept of probe particle to systems with modified inter¬ 
actions. Finally, it is a challenging and interesting prob¬ 
lem to accurately describe the residence time of ghost 
particle for any given hopping rate. 


Appendix A: Solutions for pno and pld for 
TASEP+MILK 

1. The case 5 = 7 

In the case we set 5 = 7 = I + ry, the two solutions 
to ( |39[ ) that meet the boundary conditions as derived in 
Ref. 1341 are 


Pld 


a + (1 + r]a)Q.x 
1 — (1 + r]a)r]Q.x 


(Al) 


1 - /3 + ( 77(1 - /3) + l)0(a: - 1) 
1 - ( 77(1 - /3) + l)riQ{x - I) 


(A2) 


A shock forms in the case that xi_ > x^^, where xL 
is the value for x at which pld crosses the isotherm of 
Pi = ^, and x^_^_ is the value for x at which phd crosses 
the isotherm. 


2. The case 5 = 1 + 7/;, 7 = ! — 7 /) 


The two solutions to (391 that meet the boundary con¬ 
ditions as derived in Ref. m are 


PhD = 2(1^-^) ^ for a < 1/2, 

(A3) 

^ f 2 ( 1 ^ (^0 [y{x)] + 1) + 5 for 1 - /3 > p/, 

1 2(1^ (^0 [-y{x)] + 1) + 5 for i < 1 - ^ < P/, 

(A4) 

where pld obeys the left and phd the right boundary 
condition. W[y\ is the Lambert-W function. y{x) Is given 

by 


exp 


y{x) = 
20 


1 — 7 /; 


7/ 


(2po - 1) - 1 


{x - So) + ^(2po - 1) - 1 

7/ 


, (A5) 


with po = a. So = 0 for pld and po = 1 — /3, sq = 1 
for Phd- The constant solution pMC = ^be 

equivalent of the Langmuir isotherm in the case without 
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MI. The solution pld is stable only for a < 1/2, and phd where is the valne for x at which pld crosses the 
for /? < 1/2. isotherm and is the value for x at which phd crosses 

A shock will again form in the case that x{^j^ > x\^jy, the isotherm. 


Appendix B: Linearity of R for the LD-HD phase in the large N limit 


In the large N limit, we can approximate the shock to be completely localized at Xg- The average density will then 
be equal to pld from 0 to ccg, and equal to puD from cCs to 1. We can thus approximate R as 


N 


i? =ri + n (1 


^HD 


n—2 m—2 


n—ns-\-l 


n (1 




_ m=ns + l 


(Bl) 


Here, is the on-site residence time corresponding to PLD,n-i-i and the on-site residence time corresponding to 
PHD.n-i-i- In order to evaluate this expression, we will make a few approximations. The product terms we will rewrite 
as follows: 


(1 - = 1 - wd ^ 


,LD 


771 


m=2 


m=2 


1 




\m=2 


m=2 


— CJ 


D 


3! 




.LD\ 


\m—2 


m—2 


\m—2 


-k... 


0 ^ \m—2 / j—2 ' m—2 \m—2 


J-2 


(B2) 


© 


And similarly for the prodnct over 1 — wur™. The second term in each of the sqnare brackets ensnres that the 
quantity in square brackets contains only products of terms for which m ^ m!. The factors 4, and ensure there 

is no donble counting. 


We will now investigate the components of eq. (B21 separately. The snm over m in © we can approximate as 


= V _^ 

/ y v m > / y T _ fm 


m=2 




m—2 


Jo i-i^n + ar~ n 

Furthermore, the first sum over m in © can be approximated as 

^ nn -I 

m=2 Jo " - ' - 




1 — a 


(B3) 


(B4) 


We will now look at how the different terms in eq. (|B2[) scale with N. As wu = ^ in the mesoscopic scaling, we see 


rdm = 


N 


-(^0 + a))2 n[l-{r^n + a)y 


that X © scales as 7V°, while x © scales as N Thns, in the large N limit we can discard © in eq. (B2). 
We are then left with the following expression for eq. (B2): 




m—2 


t=i 


— In » 


l-m+a) 


N 

1 — a 


1 0 


+a) 

1 — a 


(B5) 


In the second step we nsed the identity lim„_>oo J2j=i 

A similar approximation can be made for n:=2(i — oJDrJJ^), the resnlt is 




m=ns + l 




(B6) 
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We will now insert the approximations for the product terms into eq. (B11 and convert the sum over n to an integral. 
Evaluating this integral yields the final result 


R = 


Us 

1 — a 


(a - 

I i-a) + 


(B7) 


which is indeed predominantly linear in the number of sites N (note that in the mesoscopic scaling, rig is a constant 
fraction of N). An evaluation of eq. ( |B7[ ) for explicit choices for the model parameters is plotted together with the 
results from Monte Carlo simulation in Fig. 
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